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Barium titanate BaTiOs (BTO) exhibits a large electric polarization at room temperature 
and has a high piezoelectric coefficient. Therefore, it is an important compound used in various 
applications among substances with a perovskite structure. BTO undergoes three successive 
phase transitions at T ~ 403 K, 278 K and 183 K from the cubic paraelectric phase at 
high temperature to three ferroelectric phases at low temperature; the three phases exhibit 
tetragonal, orthorhombic, and rhombohedral symmetries. Many researchers have studied the 
fundamental properties and applications of BTO since its discovery. The complicated order 
process attracts interest of many theoretical physicists. However, the property of the ordered 
phases and the process of ordering of BTO are not sufficiently comprehended. To comprehend 
the ordering process, several studies were carried out to determine whether this substance 
belongs to the order-disorder or the displacive type. In a recent numerical study, it was 
suggested that BTO was not described as a simple displacive model and had order-disorder 
character.^) 

If this substance is of the order-disorder type, we expect that the ordering process mech- 
anisms describe an order-disorder model as a spin model. The order-disorder model of this 
substance was proposed for the first time by Mason and Matthias. Their model was very 
simple and was incapable of explaining the phase transitions and ordered phases of BTO. 
Subsequently, many models were proposed, but were incapable of understanding the mecha- 
nisms of the ordering process of BTO. In particular, Inoue's explanation of his experimental 
results using the Yamada model was of interest. The Yamada model quantitatively was in 
good quantitative agreement with the experimental results. ^-^ 

However, the existence of some randomness in intermediate phase is reported by an X-ray 
analysis.^-* In addition, an experimental result of hyper-Raman scattering suggests the exis- 
tence of displaced clusters at a higher temperature than the cubic-tetragonal phase transition 
temperature in the [111] directions. The previous models that include the Yamada model 
cannot describe the nature of the ordered phases. On the other hand, most of the previous 
models are used to describe only cubic-tetragonal transition. It is expected that a model 
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describing BTO over the entire temperature region will be developed. 

Therefore, the aims of the present study are to understand the phase transitions and 
ordered phases by using the order-disorder type microscopic model. We qualitatively describe 
all the ordered phases and phase transitions by our model very well. 

We consider the following Hamiltonian: 

'H = Y,^n,,ny (1) 

Here, the summation runs over the nearest-neighbor pairs of a simple cubic lattice, m 

is the spin variable that can assume eight states as shown in Fig.l and the energy parameter 
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(2) 



where Pj is defined by the density of the z-th state p{i) as follows 

Px = {p{l) + P{2) + p(5) + p{6)) - (p(3) + p(4) + p{7) + p{8)), 
Py = (P(2) + p(3) + p{6) + p(7)) - + p(4) + p{5) + p(8)), 
Pz = (P(5) + p(6) + p(7) + p{8)) - {p{l) + p{2) + p{3) + p{4)). 



(9) 
(10) 

(11) 



Pi does not represent polarization because our model is simplified, but this corresponds to 
the polarization qualitatively. In Fig. 2, we show the energy level of a neighboring spin pair 
with 6 = 0. A9 denotes a relative angle between the neighboring spins. The energy of the 
neighboring spin pair is a with = 2tan~^(l/\/2) and P with = 2tan~^(v^), namely, 
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the energy is determined by the relative angle such as in the generalized clock model/' The 
parameter 6 denotes an effect of the relative angles change by a lattice distortion. The energy 
levels modulate by 6 on the polarized phase. This model exhibits cubic symmetry as indicated 
by 5 7^ as well as 5 = 0. 

We treat this model by a pair approximation and solve these self-consistent equations 
numerically.^' In Fig. 3, we show the density of the each state at a = 0.2, (3 = 0.4, and 
S = 0.04. These values of parameters are chosen so that the ratios of the transition temper- 
atures accorded with an experiment results approximately. However, it must be noted that 
those values do not have a meaning quantitatively because our model is simplified. We observe 
three first-order phase transitions and four phases. We denote these phases as phase I, phase 
II, phase III, and phase IV from the high temperature side. Phase I is a paraelectric phase 
with cubic symmetry. Phase II, with tetragonal symmetry, has four states sharing one face 
among eight states indicating mixing. Phase III, with orthorhombic symmetry, has two states 
sharing one edge among eight states indicating mixing. In phase IV, with trigonal symmetry, 
the system completely orders one state among eight states. In Fig. 4, we show the relationship 
between the mixing states and polarization direction on each phase. These results are in good 
agreement with an experimental result. The intermediate phases correspond to the incom- 
pletely ordered phase (lOP) in the generalized clock model^) or partially disordered phase in 
a stacked triangular Ising model. ^'^' When a, which is the energy between the neighboring 
states, is less than the temperature, it is possible that these intermediate phases are realized 
for entropy gain by the mixing of the neighboring states. 

When the temperature is higher than the cubic-tetragonal transition tcmpcratme, the 
eight states are mixed and the domains of each state grow up to a specific size for the corre- 
lation with temperature nearing the transition point. Therefore, it is likely that the displaced 
clusters are detected by hyper-Raman scattering. On intermediate phases II and III, each 
of the four and two states constitute clusters of limited size; further, they are mixed. It is 
possible that the result of X-ray analysis is described by the above mentioned cluster. 

In summary, our model has three phase transitions and four phases. The symmetry and 
directions of the polarizations of the ordered phases agree with the experimental results of 
BTO. The unexplained experimental results of the intermediate phases by X-ray analysis and 
hyper-Raman scattering may be explained by our model. The intermediate phases (phases II 
and III) are considered as a type of lOP. It may be noted that the process of ordering the 
same as that was suggested by Takahashi.^^^ Although our model quantitatively describes the 
phase transitions and ordered phases, a quantitative argument is necessary for the complete 
understanding of BTO. Therefore, it may be necessary to consider the spin-phonon system 
such as in the Yamada model. Hence, there is scope for further investigation. 
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Fig. 1. States of the spin variables with a cubic symmetry. Each spin assumes eight states corre- 
sponding to the vertex of a cube. 
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Fig. 2. Energy level of a neighboring spin pair with 5 = 0. The first and second excitation states have 
a three-fold degeneration for symmetry. 




Fig. 3. Temperature dependence density of the each state at a = 0.2, /3 = 0.4, and 8 = 0.04. There 
are two intermediate phases except the paraelectric phase (phase I) and completely ordered phase 
(phase IV). On the high temperature side intermediate phase (phase II), the four states that share 
one face are dominantly well mixed. On the low temperature side intermediate phase (phase III), 
the two states that share one edge are dominantly well mixed. 




P = P//[001] P//[011] P//[111] 

(a) Phase I (b)Pliasen (c) Phase IE (d) Phase IV 



Fig. 4. The schematic figure of the relationship between the mixing states and polarization direction 
of each phase. 
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